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Three-coupled chains of quarter-filled Hubbard model with dimerization have been studied for 
periodic boundary conditions. By use of a renormalization group method, it is found that the 
interchain hopping is renormalized to zero leading to confinement when charge gap becomes 
larger than the effective interchain hopping energy. Such a result of confinement is compared 
with that of two-coupled chains. 
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Spin gap and charge gap have been studied extensively 
in coupled-chain systems. Theoretical study in terms of a 
renormalization group (RG) method shows the spin gap 
in two-coupled chains for small but relevant interchain 
hopping [Q where the electronic state is given by d x n_ y i- 
like superconducting state. The spin gap has been also 
examined for three-coupled chains [||||, ||, whose state 
depends on boundary conditions in the transverse direc- 
tion. In case of small interchain hopping at half-filling, 
there are both charge gap and spin gap for open bound- 
ary conditions (OBC) while charge gap vanishes for peri- 
odic boundary conditions (PBC). A system with infinite 
chains at half-filling exhibits charge gap but no spin gap 
1- 

Although these works treat the relevant interchain hop- 
ping, a theory also suggests the irrelevant one-particle 
interchain hopping for half- filling ||. In order to un- 
derstand optical experiments on quasi-onc-dimensional 
organic conductors ||[7|] indicating a confinement - de- 
confinement transition, two-coupled chains at half-filling 
have been studied where a transition from deconfinc- 
ment to confinement occurs for umklapp scattering being 
larger than a critical value ^,^]. In the present paper, 
such a transition is examined for three-coupled chains. 
Although the results of two chains seem to explain the 
experiments, it is not yet known if such a transition can 
be expected for many chains. 

The Hamiltonian for quarter-filled three-coupled chains 
with on-site repulsion (U) and interchain hopping (t±) 
is written as 

n = - E V + (ck, c j+lal + h. c 

-t±^2 (c\ al c 3al+1 + h.c.) + n m n 3iii l ) 

3,1 



3,°,l 



where rij^i — c^ al Cj a i and Cj a \ denotes the annihilation 
operator of the electron at the j-th site of the Z-th chain 
(I =1, 2, 3) with spin tr(=Ti4)i an d Cjai — £jo\- In 
Eq. (1), t d denotes dimerization along the chains and 
the case for PBC is studied. 



Diagonalizing the id-term, which leads effectively to 
the half-filled band Q , we consider an effective Hamilto- 
nian H d consisting of the lower band with fermion opera- 
tors of dka-i ■ The terms for the interchain hopping can be 
diagonalized by introducing Fourier transform, a^ap, = 
(I/n/3) J2i=i ex P MM/^) (I dkd with ky(ji) = (27r/3)/i 
(/i = 0, ±1). The kinetic term is written as H^ = £ fc a u 



a 
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with e(k,k y 



t 2 cos 2 ka + t\ sin 2 ka- 



cosfc,,, which is rewritten, in terms of the linearized 



dispersion, as H d K = £ fciP>CT|A1 v F (pk - k Ffl ) a kpt7fl 
with p being the index of the branch p = + (— ) corre- 
sponding to right moving (left moving) electrons. Fermi 
momenta are given by k F o = k F + 2t±/v F and k F ± — 
k F - t±/v F where v F = V2ta [1 - (t d /t) 2 } j^J\ + (t d /t) 2 
|J. Following the conventional g-ology, coupling con- 
stants are given by g 1 ± = g 2 ± = = Ua, g$ = 
Ua(2t d /t)/[l + (t d /t) 2 ] § and 9lll = g 2ll = g 4ll = 0. 

Applying the bosonization method, we introduce phase 
variables 9 PP and 9 ap expressing fluctuations of the charge 
density and spin density for the /i-band |[o|], where the 
conjugate phase is introduced by 4> v i ^(x 1 )]- = 

m5 v y S^^i sgn(a; — x') (y,v' — p,cr). Thus we obtain 
the total Hamiltonian given by H d = H + Hi where 
Hq expresses bilinear terms of density operators, and Hi 
denotes nonlinear terms. Here we define the new phase 
variables as Y x = (X 1 +X 2 +X 3 )/V3, Y 2 = {X x -X z )/>/2 
and Y 3 = (-Ax + 2A 2 - A 3 )/V6 where (Yi,F 2 ,r 3 ) = 
(0 P ,0ci,0c2) and (Xi,X 2 , A 3 ) = (<9 p+ , 9 pQ , 6> p _ ) and the 
same transformation is applied to p — > a and also to 
the conjugate variables cj> vll . Thus Ho is rewritten as 
Wo = £„(W47r) / dx[K- x (d6 v f+K v {d<t> v f\ with v = 
p,Cl,C2,cr,Sl,S2, where v p(u) = v F [l + {-)U/ttv f } 1 / 2 , 



K, 



Ua/nv F ] 1/2 , u i = v S i = v C 2 = v S2 = 
v F and Kqi — ATsi = Kq2 — K§ 2 = 1. The Hamilto- 
nian, Hi is divided as Hi = Hi + Hz + H 2 + Hn where 
the respective term is written as (i — 1, 3, 2, ||) 



Hi = — ^ Giz / dx cos ®iz ■ hiz 



(2) 



(Z = A ~ Z) and a is a cutoff of the order of lat- 
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Fig. 1. Quantities tj_(i)/t, K p (l) and K a (l) as a function of I for Fig. 2. Normalized quantity tj_(l)/t for choices of U/t = 3, U c /t 

U/t = 4, t±/t = 0.1 and t d /t = 0.1. The inset shows tj_(l)/t, and 4.5 with t x /t = 0.1 and t d /t = 0.1 where t/ c /t ~ 3.7. The 

K p (l), Kqi(2), Kq2(3), K a (&), Kgi(5) and ii"s2(6) for U/t = 4, solid and dotted curves are explained in the text. 
t ± /t = 0.1 and t d /t = 0.025. 



tice constant and hiz denotes the product of the Majo- 
rana fermion operators introduced to retain the anticom- 
mutation relation of the field operators H . The phase 
variables of the backward scattering term with opposite 
spins, Hi, are given by Q 1Z = (2/V3)6 a + e lz where 
Q 1A = eV20si - Oib = 2^/38 S2 , Q lG = 

[(>'0si +0 S2 /V3) + e{6 C i - e'V3(0 C 2 + 2V6t A _x/v F ))]/V2, 
Qw = eV29 C i - ^2/3^2, Qie = [(e'^si + fa/Vty + 
e(0 sl -eV30s 2 )]/\/2, B 1F = ev^si- v^sa, § 1G = 
(^2+e , V30ci+e30c2-ee , V3</>si)/v / 6ande m = [(e"0 si 
0s 2 /\/3) + e' (0 C i + e"V3(9 C 2 + 2<ft>t x x/v F )) +e(e"30 C i - 
v ? 3(/'C2)-ee'(0si+e"\/30s2)]/2\/2. In Eq. (2), the sum is 
taken implicitly with respect to e, e' and e"(= ±) which 
lead to the distinction for hiz but not for Giz- The 
umklapp scattering terms, TI3, are obtained from 7i\ by 
replacing G 1Z -> G 3Z and (9 a , 9 S i, S2 ) «-> (6> p , 6» C i, #02 + 
2v6t_Lx/up). The phase variables of the forward scatter- 
ing term with opposite spins, H 2 , are expressed as Q 2 e = 
[(0ci ~ e' V3(9c2 + 2V6t ± x/v F )) + e{<j> S i - e' V3^ S2 )]/ A 
9 2 f = V20ci+ev^^si, e 2G = (e C i-e'V30s2-ee'y3</>c 2 
e^ S i)/^and Q 2H = [{9 cl +e"V3(6 C 2 + 2V6t ± x/v F ))- 

e'(e"30si-v / 30 S2 )-ee'(e"30ci-\/3^c 2 )+e(0si+e"v / 3</'s 2 
The forward scattering term with parallel spins, Tin is 
obtained from H 2 by replacing G 2 z — > G\\z(Z = E ~ 
iJ) and (#si,#s 2 ) (0Si,0S2)- Coupling constants are 
given by GiA = G iB = G,g = G,.d = G l£; = G iF = 
Gjg = G iff = gi/6irv F (i = 1,3) and Gj B = G,i? = 
G,g = G iff = 3i/67ri; F (t = 2, ||) where = g 2 \\ - g x \\. 
Non-linear terms of forward scattering with the same p 
branch are discarded because these effect is negligibly 
small. By assuming scaling invariance with respect to 
a — ► a' = ae dl , the second order RG equation for the 
interchain hopping is given by 



where t± 



~ t± = h - F (t x , {G tZ }) K C2 , 
= ^/(.p^ 1 ), F(t ± ,{G lZ }) = [G\ c 



(3) 
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GL]Ji(6fx) + [G? ff +G! ff +Gjj ff ]J 1 (3fx) + [Gi 
Gl F ]J 1 (4tx)/3+[G| c +G| E +G| G ]J 1 (2?x)/3+G| B J 1 (8fx)/3+ 
Gg ff Ji(i^)/3 and J„ is n-th Bessel function. Equa- 
tion (3) is solved together with RG equations for Giz 
where I is related to energy scale ui or temperature T 
by I = ln(W/w) or ln(W/T) with W(= w F a _1 ) being of 
the order of band width. We take a = 2a /ir ||. It is 
noted that the r.h.s. of eq. (3) for small t± is reduced to 
t±[l - (G\ + G\ + G| + G|)/2] which becomes the same 
as that of many chains [^| . 

In Fig. 1, the /-dependence of coupling constants is 
shown for t d /t = 0.1 and t d /t = 0.025 with U/t = 4 
and t±/t = 0.1. For t d /t = 0.025 (inset), the inter- 
chain hopping t±(l) becomes relevant corresponding to 
dcconfincment where K p remains finite and Kqi, K a 
and Ks 2 (Kq2 and Ksi) decrease to zero (become in- 
finite). The curves are calculated, for simplicity, by set- 
ting J n (l) = for I > l n with l n corresponding the 
first node of the Bessel function, although such a treat- 
)]/a^It gives negligible difference in the numerical results. 
When K v (l) decreases to zero or increases to infinity, the 
corresponding phase is locked leading to a formation of 
gap, where relevant coupling constants are Gib(— > —00), 
Gid(-> -00), Gif(— > -00), Gig(— > +00), G 2F {-^ -00) 
and G 2 g{-^ +00). There are two kinds of gap for charge 
fluctuations and three kinds of gap for spin fluctuations. 
For td/t — 0.1 (main figure), one finds confinement where 
t±(l) decreases to zero after taking a maximum. The 
/-dependence of K p (l) implies charge gap in the total 
charge fluctuation, which is in contrast to the case of de- 
confincmcnt. For t d /t = 0.1, exponents K v (l) (y =C1, 
C2, (J, SI, S2) remain finite, i.e., K v (l) ~ 1 at I ~ 3. 
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Fig. 3. The critical values U c as a function of tj_ with fixed 
td/i = 0.05(solid curves) and 0.1(dashed curves) for three- 
coupled chains and two-coupled chains |j| . The inset denotes the 
log t<f /t_L-log tj_ plot for t d /t = 0.0 and 0.1 with fixed U/t = 4.0 
where the arrow shows the critical value for the confinement. 



In Fig. 2,t±(l)/t is shown with the fixed U/t = 3, U c /t(~ 
3.7) and 4.5 where solid (dotted) curves denote t±(l)/t 
for < I < / A (I > I a) with K p (l A ) = K p /2, The case 
for I > /a is invalid since the magnitude of the umk- 
lapp scattering increases to infinity. In the curve t±(l), 
there is a maximum given by max[t±(l)/t] ~ 3.4 (0.5) 
at I ~ 4.0 (2.2) for U/t = 3.0 (4.5) although the max- 
imum for U/t — 3.0 is located in the region of I > 
Thus the case for U/t = 3.0 (4.5) corresponds to decon- 
finement (confinement). The boundary between confine- 
ment and deconfincmcnt is determined by the condition 
that max.[t±(l)/t] ~ 1 at U = U c . 

In Fig. 3, the independence of U c is shown for td/t = 
0.05(solid curves) and 0.1(dashed curves) where confine- 
ment (deconfinement) is obtained for U > U c (U < 
U c ). For comparison, the corresponding results for two- 
coupled chains ^ are also shown. The critical values 
for three-coupled chains is smaller than that for two- 
coupled chains. Since the RG equation of two-coupled 
chains corresponding to Eq. (3) includes the Bessel func- 
tion with only J\(8i±), the effect of umklapp scattering 
for three-coupled chains is stronger than that for two- 
coupled chains. The effective interchain hopping tj° can 
be evaluated from t°f = iexp[— l c g] where t±(l c g)/t = 1 
pi]] . In the inset, tf± /t± is shown as a function of 
t± on a logarithmic scale. The power-law behavior of 
t c f'° (= ^(53 = 0)) for small t± is consistent with 
the analytical formula tf '° ~ t ± (t J _/W) a °/ < - 1 - a ^ with 
a = {K p + K- 1 + K a + K- 1 - 4)/A §. In the pres- 
ence of the dimerization, t c f is reduced from that of the 
power-law behavior and has a critical value below which 
the confinement occurs. 

The confincmcnt-dcconfincmcnt transition is examined 




Fig. 4. The phase diagram of confinement (I) and deconfincmcnt 
(II) in case of td/t = 0.05 (solid line) and 0.1 (dashed line) on 
a plane of the effective interchain hopping t^'° and A^ D . The 
inset denotes the corresponding phase diagram on the plane of 
t± and AJ D . 



in terms of the charge gap induced by the umklapp scat- 
tering. The charge gap for the single chain is obtained as 
the function of U and id by a method A* = W exp[— I a] 
with K p {l&) = K p /2 for t± = 0. Such gap has a mean- 
ing of a characteristic energy of the umklapp scatter- 
ing even for the deconfincd region in which the charge 
gap is reduced to zero due to the presence of the mis- 
fit for all the nonlinear terms of umklapp scattering. 
In the inset of Fig. 4, a phase diagram of confinement 
(I) and deconfinement (II) is shown on the plane of the 
bare interchain hopping t± and Ap D where the bound- 
ary for two-coupled chains || are also shown for com- 
parison. Although the ratio of A* D /t^ is nearly con- 
stant, the curve is slightly convex upward for small t±. 
In the main figure of Fig. 4, the phase diagram with 
the same parameter is shown on the plane of t e ± '° /t and 
A^ D /t. The quantity t c f'° denotes the effective inter- 
chain hopping, which is renormalized by the intrachain 
interaction without umklapp scattering. The ratio of 
Ap B to t c f'° at the boundary is estimated as follows 
when 0.05 < t±/t < 0.3. The ratio for three-coupled 
chains is given by 1.0 ^ A* /^ 1 ' ^$ 1.1 while that for 
two-coupled chains is given by 1.8 ^ A^ D /t c f'° ^ 1.9. By 
noting that the id-dependence of the boundary is very 
small, it turns out that the confinement-deconfinement 
transition is determined by the competition between the 
charge gap and the effective interchain hopping energy. 

We briefly discuss the case for OBC where the RG 
equation for the interchain hopping takes more compli- 
cated form due to twelve coupling constants for umklapp 
scattering. The small difference between PBC and OBC 
is expected since the RG equation for the interchain hop- 
ping has the same limiting form as that of PBC for small 
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i± . Actually, we find almost the same boundary as Fig. 3 
when the solution of single chain is substituted for the 
RG equation of interchain hopping. 

In summary, the confinement-deconfinement transition 
in the three-coupled Hubbard chains with dimerization 
has been shown for PBC when the effective interchain 
hopping energy becomes of the order the charge gap in- 
duced by the umklapp scattering. 
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